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Abstract. We have made new abundance determinations 
for a sample of NaMgAl stars. These stars, which are 
a subgroup of the nearby metal-rich field F and G disk 
dwarfs, were first identified by Edvardsson et al. (1993) 
on the basis of their apparent enrichment in Na, Mg and 
Al relative to other elements. The discovery of a plan- 
etary companion to the nearby solar type star 51 Peg 
(Mayor & Queloz 1995) combined with Edvardsson et al.'s 
earlier identification of 51 Peg as a NaMgAl star high- 
lighted the group's potential importance. Our new analy- 
sis, which uses new spectra of higher resolution and bet- 
ter wavelength coverage than the analysis of Edvards- 
son et al., shows that the Na, Mg and Al abundances 
of the NaMgAl stars are indistinguishable from those of 
non-NaMgAl stars with otherwise similar properties. The 
group thus appears to be spurious. 

Our study, which includes 51 Peg, also provides the 
most complete set of abundances for this star available to 
date. The new Fe abundance, [Fe/H] = +0.20 ± 0.07, of 
51 Peg confirms earlier measurements of its metal richness. 
Abundances for 19 other elements, including C, N and O, 
reveal a fairly uniform enrichment similar to that of Fe 
and show no evidence of abnormality compared to other 
metal rich stars of similar spectral type. 

Key words: Stars: abundances ~ Stars: fundamental pa- 
rameters - Galaxy: abundances ~ Galaxy: evolution 



1. Introduction 

In a survey of abundances in nearby field F and G dwarfs 
Edvardsson et al. (1993, EAGLNT below) suggested the 
existence of a new group of stars that are enriched in Na, 
Mg and Al relative to other elements. EAGLNT identi- 
fied eight of these "NaMgAl" stars in their survey of 189 



stars. The average enhancements relative to other stars of 
similar metallicity are AIX/Fe] = 0.12,0.14,0.10 for Na, 
Mg and Al, respectively. El These stars are also character- 
ized by metal richness - all have [Fe/H] > 0.05. EAGLNT 
suggested that perhaps "the stars have been affected by 
mass transfer from a now 'dead' companion" . The recent 
discovery of planetary companions to nearby solar type 
stars has given this speculation unexpected immediacy. 

The first of these discoveries was 51 Peg (Mayor & 
Queloz 1995). Mayor and Queloz's radial velocity obser- 
vations revealed the presence of a companion with a min- 
imum mass A^sini = 0.47 ± 0.027Wj (A^j — the mass 
of Jupiter), a result confirmed by Marcy et al. (1997). 
A remarkable feature of the companion is its short or- 
bital period - only 4.2 days. This places the companion 
only 0.05 au from its parent star, which is puzzlingly close 
for a planet of Jupiter- type mass. 51 Peg was included in 
EAGLNT's survey and labelled as a NaMgAl star. The 
suggestion as to these stars' origin thus raises the possi- 
bility that 51 Peg's companion may be an exotic form of 
stellar remnant, rather than a planet. 

Further discoveries of planetary companions have fol- 
lowed 51 Peg. Three of these additional new systems - 
V And, T Boo and Cue (Butler et al. 1997) - mimic 
51 Peg in that they are solar type stars and their com- 
panions are very close to their parent stars. Two of these 
systems (r Boo and 55 Cnc) were not included in the sur- 
vey of EAGLNT, while the third {v And) was, but was not 
identified as a NaMgAl star. This extremely limited sam- 
ple thus suggests that the possible correlation between the 
NaMgAl stars and close planetary companions is imper- 
fect, nonetheless it is clear that the NaMgAl stars deserve 
closer examination. 

Here we report an abundance analysis of the NaMgAl 
stars based on new observations made with the McDon- 
ald Observatory's 2.7m telescope and 2dcoude spectrom- 
eter. The generous wavelength coverage of this instrument 
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Table 1. Programme star data. Column 3) indicates whether the star was identified as a NaMgAl star in EAGLNT. Columns 
4)-7) show the model atmosphere parameters, 8)-10) show velocities relative to a local standard of rest where the U velocity is 
taken positive towards the galactic centre. Columns 11)-15) give kinematic data and ages adopted from EAGLNT 
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Fig. 1. Two All and two Fei lines near 6700 A in HR448, a NaMgAl star, and \x Cnc, a non-NaMgAl star. The spectrum of ji 
Cnc has been wavelength shifted to match that of HR448. The two stars have similar effective temperatures and gravities (see 
Table 1) so the striking similarity of their spectra indicates that any enhancement of Al relative to Fe in HR448 is very small 



provides several more lines of the three elements than 
EAGLNT were able to use; their observations only cov- 
ered a few limited wavelength regions. And, for the north- 
ern hemisphere stars in the EAGLNT survey, the high 
resolution {R =60000) of this instrument also provides 
more accurate equivalent widths. The high resolution is 
particularly significant for analysis of metal-rich stars like 
51 Peg. (The northern hemisphere stars of the survey of 
EAGLNT, were observed with a resolution of 30 000, while 
those of the southern hemisphere were observed with reso- 
lutions of 60 000 or 80 000). The wide wavelength coverage 
also allows access to some additional elements not included 
previously. 



The main objective of the new study of the NaMgAl 
stars, therefore, is to inspect the enhancements of Na, Mg 
and Al, which are small, more closely and to see if any 
other elements are involved. A secondary objective is to 
expand the set of elements with measured abundances for 
51 Peg. 

2. Observations 

The observations were made at McDonald Observatory 
with the 2.7 m telescope and 2dcoude echelle spectrometer 
(Tull et al. 1995). The detector was a Tektronix CCD with 
24 /xm^ pixels arranged in a 2048x2048 pixel format. 

We observed six of the eight NaMgAl stars identified 
by EAGLNT and three other stars in the EAGLNT sur- 
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vey not identified as NaMgAl stars, but with similar effec- 
tive temperatures, gravities and metallicities. These nine 
stars, which include 51 Peg, are given in Table 1. Finally, 
we observed an asteriod (Iris) in order to provide a solar 
spectrum recorded under the same circumstances as the 
stellar spectra. 

We observed the wavelength interval 4000 - 9000 A 
approximately. The coverage is complete from the start 
of this interval to 5600 A and substantial, but incomplete, 
from 5600 A to the end of the interval because of gaps 
between the end of one echelle order and the beginning 
of the next. The slit width was set to project onto two 
pixels, which gave a resolution of 60 000. From ^ 5500 
to '-^ 9000 A the extracted one-dimensional stellar spec- 
tra have a typical signal-to-noise ratio of ~ 400, while at 
shorter wavelengths than ~ 5500 A the signal-to-noise ra- 
tio decreases with decreasing wavelength because of the 
decline of the stellar (and flat field) fluxes towards shorter 
wavelengths. 

Figure 1 shows the two Al i lines near 6700 A in a Na- 
MgAl star and a non-NaMgAl star. The two stars have 
similar effective temperatures and gravities so differences 
in the strengths of their lines primarily reflect abundance 
differences. The strengths of the Al i lines relative to the 
Fe I lines in the two stars are very similar suggesting that 
any enhancement of Al relative to Fe in the NaMgAl star 
is small. 

The data were processed and wavelength calibrated 
in a conventional manner with the IRAF package of pro- 
grams on a SPARC 5 workstation in the Astronomy De- 
partment at the University of Texas. Lines suitable for 
measurement were chosen for clean profiles, as judged by 
inspection of the solar spectrum at high resolution and 
signal-to- noise ratio (Kurucz et al. 1984), that provide re- 
liable equivalent widths in all, or most, of the programme 
stars. Moore et al. (1966) was our primary source of line 
identification. The equivalent width of each line was mea- 
sured with the IRAF measurement option most suited to 
the situation of the line; usually this was the fitting of a 
single, or multiple, Gaussian profile to the line profile. 

Table 2 gives the list of lines measured. The list in- 
cludes most of the lines used by EAGLNT and many ad- 
ditional ones. In particular our new list contains more Na, 
Mg and Al lines than theirs; specifically three Nai, four 
Mg I and seven Al i lines in the new line list compared with 
two Nai, two Mgi and two All lines in theirs. The new 
line list also has more lines of most other elements and it 
includes some new elements, notably C and N which are 
represented by six C i lines and two N i lines. 

Figure 2 compares our new equivalent widths with 
those of EAGLNT. The overall agreement is good. In a 
few cases, however, we note a tendency for the new equiv- 
alent widths to be slightly larger than those of EAGLNT; 
51 Peg is the most marked case. This can be attributed 
to the fact that for most of the stars in the programme 
the resolution (60 000) of the new observations is higher 



than that (30 000) of EAGLNT's observations thus leai 
ing to a better defined and slightly higher continuum. □ 
The behaviour of 51 Peg, which as one of the coolest and 
most metal rich of the stars has one of the most crowded 
spectra, supports this explanation. 

3. Analysis 

The LTE abundance analysis was made strictly relative 
to the Sun and adopts the same procedures, software and 
model atmospheres as were used by EAGLNT. 

3.1. Line data 

Basic line data for all the measured transitions were ob- 
tained from the VALD atomic line data base (Piskunov et 
al. 1995) by the "show line" command. Astrophysical oscil- 
lator strengths were determined for each line by requiring 
that the equivalent width derived from the solar model 
be equal to that measured from the solar fiux spectrum 
as reflected by Iris. The parameters of the solar model 
are given in Table 1 and the solar chemical abundances 
were adopted from Grevesse et al. (1996). New quantum 
mechanical calculations of pressure line-broadening prop- 
erties by Anstee & O'Mara (1995) and Barklem & O'Mara 
(1997) were adopted for several lines of s-p, p-s, p-d and d- 
p transitions in neutral species. There is a small systematic 
difference between our solar equivalent widths and those 
measured by EAGLNT. The weakest lines are in the mean 
about 2% stronger in EAGLNT, and the strongest lines 
about 2% weaker. This difference/uncertainty (together 
with the line-broadening mentioned above) is reflected in 
the gf values. However, since the programme stars are ob- 
served and measured with identical procedures, we neglect 
the very small differential errors which may be caused by 
these uncertainties. 

A few lines were excluded from the analysis after a 
first abundance analysis of all stars, since they give abun- 
dances that deviate strongly, often as a function of ef- 
fective temperature, from those of other lines of the same 
species: Sci 6239.36, Til 5299.98, Vi 6150.15, Fei 5127.37, 
Nil 5099.94 and 7788.15 and Yi 6687.51 A. For some of 
these, asymmetries could be traced in the spectra, while 
for others we suspect undetected blends or line misidenti- 
fications. For the 7788.15 A Ni i line a laboratory oscillator 
strength is quoted in the VALD data base which is 0.5 dex 
lower than our solar value, which means that most of the 
observed line absorption in the Sun may be caused by 
another transition. Therefore, these lines were not used. 
Also the Sen 5318.36 A line was found to be asymmetric 

^ Of the nine stars in our programme seven {fi Cnc, 20 LMi, 
24 LMi, 9 Com, 11 Aqr, HR 8472, 51 Peg) had been observed 
from McDonald Observatory at a resolution of 30 000 and two 
(HR 448, HR 1536) had been observed from ESO at a resolu- 
tion of 80 000 in EAGLNT's study. 
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Table 2. Atomic line data and solar equivalent widths. The 
columns give: 1) the line rest wavelength, 2) the excitation en- 
ergy of the lower level of the transition, 3) the logarithm of the 
product of the statistical weight and the oscillator strength, 4) 
an enhancement factor to the classical van der Waals damping 
constant; An "A" or a "B" indicates that the pressure damping 
was instead calculated according to Anstee & O'Mara (1995) or 
Barklem & O'Mara (1997) respectively, 5) the radiation damp- 
ing constant, and 6) the equivalent width measured from the 
solar spectrum as reflected by the asteroid Iris 



Table 2. Continued 



A 


Xi 


log gf 




Trad 


Wa0 


[A] 


[eV] 






s-i 


[mA] 


C I loge© = 


8.55 








5380.32 


7.69 


-1.68 


2.5 


1.710® 


22.3 


6587.62 


8.54 


-1.14 


2.5 


3.910® 


14.6 


7111.45 


8.64 


-1.18 


2.5 


1.810® 


11.2 


7113.17 


8.65 


-0.78 


B 


2.5 10® 


22.7 


7115.17 


8.64 


-0.68 


B 


4.0 10® 


26.8 


7116.96 


8.16 


-1.30 


2.5 


1.410® 


18.5 


N I 


log £0 = 


7.97 








7468.27 


10.34 


-0.06 


A 


3.110^ 


4.1 


8629.16 


10.67 


0.08 


A 


3.0 lO'^ 


3.4 


O I 1 


log £0 = 


8.87 








6158.17 


10.74 


-0.37 


B 


5.410^ 


4.2 


7771.95 


9.14 


0.36 


A 


4.810^ 


70.6 


7774.18 


9.14 


0.21 


A 


4.7 lO'' 


60.9 


7775.40 


9.14 


-0.02 


A 


4.5 lO'' 


47.9 


Na I 


log £0 = 


= 6.33 








4497.68 


2.10 


-1.59 


2.0 


6.410^ 


39.3 


6154.23 


2.10 


-1.59 


2.0 


6.410^ 


39.6 


6160.75 


2.10 


-1.32 


2.0 


6.710'' 


57.7 


Mg I 


log £0 


= 7.58 








6318.71 


5.11 


-2.08 


2.5 


2.810^ 


39.8 


6319.24 


5.11 


-2.28 


2.5 


3.0 lO-""' 


27.8 


6965.41 


5.75 


-1.75 


2.5 


3.6 10^ 


24.1 


7759.37 


5.93 


-1.68 


2.5 


3.310^ 


21.0 


Al I 


log £0 = 


: 6.47 








6696.03 


3.14 


-1.59 


2.5 


3.010* 


38.1 


6698.67 


3.14 


-1.91 


2.5 


3.010* 


21.9 


7835.32 


4.02 


-0.72 


2.5 


7.9 10^ 


48.1 


7836.13 


4.02 


-0.57 


2.5 


7.9 10'' 


61.0 


8772.88 


4.02 


-0.46 


2.5 


8.310'^ 


75.2 


8773.91 


4.02 


-0.29 


2.5 


8.3 10^ 


94.4 


8828.87 


4.09 


— 1.89 


2.5 


4.2 10® 


4.0 


Si I 


log £0 = 


7.55 








6125.03 


5.61 


-1.53 


1.3 


1.110® 


33.4 


6142.49 


5.62 


-1.47 


1.3 


8.810^ 


36.5 


6145.02 


5.61 


-1.42 


1.3 


9.8 10^ 


39.9 


6155.14 


5.62 


-0.75 


1.3 


3.6 10® 


88.7 


6848.57 


5.86 


-1.65 


1.3 


1.1 10® 


18.3 


7455.39 


5.96 


-2.01 


1.3 


4.0 10^ 


7.5 


7800.00 


6.18 


-0.73 


1.3 


3.3 10® 


55.7 


SI log 60 = 


7.33 








6046.02 


7.87 


-0.41 


2.5 


1.0 10^ 


17.3 


6743.58 


7.87 


-0.64 


2.5 


1.1 10^ 


10.3 


7686.13 


7.87 


-1.15 


2.5 


1.710® 


3.4 


K I 


log £0 = 


5.12 








5801.75 


1.62 


-1.51 


2.5 


3.1 10® 


2.0 


7698.98 


0.00 


0.01 


A 


5.8 10^ 


159.5 



A 


Xi 


logs/ 




Trad 


Wa0 


[A] 


[eV] 






s-i 


[mA] 


Ca I 


log £© = 


= 6.36 








5867.57 


2.93 


—1.64 


1.8 


2.6 10 


25.3 


6166.44 


2.52 


— 1.22 


B 


1.9 lO'' 


71.7 


6169.04 


2.52 


—0.88 


B 


2.0 lO'^ 


96.0 


6455.61 


2.52 


— 1.43 


A 


4.6 lO'^ 


58.6 


6464.68 


2.52 


—2.41 


A 


4.6 lO'^ 


12.6 


6572.80 


0.00 


—4.33 


A 


2.6 10^ 


33.6 


So I 


log £0 = 


= 3.17 








5520.51 


1.87 


0.62 


A 


9.3 10 


9.0 


So II 


log £0 


= 3.17 








5318.36 


1.36 


—1.81 


2.5 


1.4 10 


13.5 


6245.62 


1.51 


—1.10 


2.5 


4.3 10* 


36.7 


6604.60 


1.36 


— 1.25 


2.5 


1.5 10* 


37.0 


Ti I 


log£0 = 


= 5.02 








5113.45 


1.44 


—0.91 


A 


2.2 10 


27.5 


5219.71 


0.02 


—2.29 


A 


6.710® 


28.7 


5426.26 


0.02 


—3.11 


A 

A 


1.7 10® 


6.2 


5766.33 


3.29 


0.26 


B 


1.4 10* 


9.4 


5866.46 


1.07 


—0.88 


A 


6.4 10*^ 


48.7 


6091.18 


2.27 


—0.47 


A 


8.5 lO'' 


15.9 


6126.22 


1.07 


— 1.44 


A 


9.9 10® 


22.8 


6258.11 


1.44 


—0.47 


A 


1.7 10* 


52.2 


V I log £0 = 


4.00 








5727.06 


1.08 


—0.03 


A 


7.1 10 


39.4 


6039.74 


1.06 


—0.73 


A 


4.0 10^ 


13.3 


6090.22 


1.08 


-0.15 


A 


4.0 lO'^ 


34.2 


6111.65 


1.04 


—0.80 


A 

A 


3.9 lO'^ 


12.1 


6216.36 


0.28 


-0.87 


A 


3.1 10® 


38.4 


6224.51 


0.29 


—1.91 


A 


1.2 10® 


5.7 


6251.83 


0.28 


—1.42 


A 


3.1 10® 


15.8 


Cr I 


log £0 = 


= 5.67 








6330.10 


0.94 


—2.94 


A 


2.4 10 


27.1 


Cr II 


log £0 


= 5.67 








5305.87 


3.83 


—2.06 


2.5 


2.6 10 


26.6 


Fe I 


log £0 = 


= 7.50 








5067.16 


4.22 


-0.96 


B 


2.1 10 


73.5 


5090.78 


4.26 


—0.61 


B 


2.1 10* 


96.5 


5109.66 


4.30 


—0.82 


B 


2.1 10* 


79.4 


5141.75 


2.42 


—2.20 


A 


1.5 10* 


89.6 


5358.12 


3.30 


—3.20 


A 


2.0 10* 


10.3 


5809.22 


3.88 


— 1.70 


A 


5.1 lO'^ 


50.6 


5849.69 


3.69 


—2.98 


A 


5.5 lO'' 


7.8 


5852.23 


4.55 


—1.23 


B 


1.9 10* 


41.2 


5855.09 


4.61 


— 1.57 


B 


1.9 10* 


22.4 


5856.10 


4.29 


-1.60 


1.4 


8.6 lO'' 


33.4 


5858.79 


4.22 


—2.20 


B 


2.7 10* 


13.8 


5859.60 


4.55 


—0.71 


B 


1.9 10* 


72.8 


5861.11 


4.28 


—2.36 


B 


2.1 10* 


8.9 


5862.37 


4.55 


-0.49 


B 


1.9 10* 


88.4 


6151.62 


2.18 


-3.30 


A 


1.6 10* 


50.9 


6157.73 


4.08 


-1.22 


1.4 


5.0 10^ 


63.4 


6159.38 


4.61 


-1.88 


B 


1.910* 


12.6 


6165.36 


4.14 


-1.50 


1.4 


8.8 10^ 


46.3 


6173.34 


2.22 


-2.88 


A 


1.710* 


69.7 


6200.32 


2.61 


-2.39 


A 


1.010* 


75.3 


6436.41 


4.19 


-2.39 


1.4 


3.0 10'' 


10.2 



J. Tomkin et al 



Table 2. Continued 



A 


Xi 


logs'/ 




Trad 


Was 


[A] 


[cV] 








[mA] 


Fe I cont. log 


e© = 7.50 






ooyi. oo 


4.oy 


1 QQ 

— i.yy 


1 A 
1.4 


1.4 10 


10.8 


DDUO.U4 


Z.Zo 


Q QC^ 

— o.yo 


A 


1.7 10* 


18.7 


fifiQQ 1 A 
DDyy. 14 


A t^Q 

4.oy 


10 

— z.iz 


1 A 
1.4 


1.4 10* 


8.4 




'4:.OU 


— 1.44 




2.4 10* 


21.5 


D / ZO.OO 


A ^(^ 

4.iU 


01 
— Z.Zl 


A 


2.1 10* 


17.9 


D / OZ.U 1 


A ciS 
4.0o 


01 
— Z.Zl 


1 A 
1.4 


6.0 lO'^ 


7.1 


D / OO. 10 


A «A 
4.04 


1 /I 
— 1.40 


1 A 
1.4 


2.3 10* 


27.7 


DoO / .ZO 


A ns 

4.Uo 


07 

— z.u / 


1 A 
1.4 


2.5 10^ 


23.5 


/ / 01. IZ 


/I QQ 

4.yy 


SO 
— U.OU 


13 


6.4 10* 


46.5 


/ oUZ.Dl 


OQ 
o.uy 


1 

— l-OD 


£> 


6.3 10* 


16.1 


1 044.0 1 


4.o4 


1 71 
— 1. / 1 


1 A 
1.4 


2.3 10* 


12.8 


O / 4 / .44 


09 
O.UZ 


— 0.00 


A 


8.0 10^ 


18.1 


Q7C7 on 
o / / .ZU 


Z.o4 


01 
— Z.Ui 


A 
l\ 


7.7 10'^ 


97.4 


■Pc TT 




- 7 CiO 








oiuu.oo 


9 SI 
Z.Ol 


4 1 
— 4.10 


9 

Z.O 


3.4 10 


22.1 


ozoD.yo 


QQ 

z.oy 


A OS 




Z.O 


3.4 10* 


22.1 


04Z0.ZD 


Q 00 
O.ZU 


— 0.01 


K 
Z.O 


3.0 10* 


41.8 


^497 SO 


(S 79 

U. 1 z 


-L .UO 


Zi.O 


3.5 10* 


4.4 


D14y .zo 


Q SQ 

o.oy 


7S 
— Z. / o 




Z.O 


3.4 10* 


37.0 


CO/17 t^^i 
DZ4 / .00 


o.oy 


0/11 
— Z.41 


K 
Z.O 


3.4 10* 


54.4 


OOOy .40 


9 SO 

Z.oy 


A 1 
— 4.10 


9 K 
Z.O 


2.9 10* 


19.9 


OOoO. / Z 


0.00 


Oil 


9 

Z.O 


4.1 10 


9.8 


dA QO 
D4oZ.Oo 


9 SO 

Z.oy 


Q A/1 
— O.D4 


9 

Z.O 


2.9 10* 


41.6 


D40D.oy 


^ QO 

o.yu 


90 
— Z.ZU 


9 

Z.O 


3.4 10* 


64.7 


7/1 70 70 

/ 4 / y. ( U 


Q so 
o.oy 


Q (^7 
— 0.0 ( 


9 Pi 
Z.O 


3.1 10* 


9.3 


i 010. o4 


Q 00 

o.yu 


f^O 
— O.OU 


9 K 
Z.O 


4.1 10* 


12.6 


78/1 1 '^7 


00 

o.yu 


A 04 
— 4.U4 


9 

Z.O 


3.0 10* 


4.3 


tm; t 

IM 1 1 




o.zo 








OUoZ.oO 


O.DD 


— u.oy 


R 
Jd 


1.8 10 


67.9 


0Uo4. 1 1 


(^S 
O.Do 


1*^ 
— U. 10 


"R 
13 


1.3 10* 


94.8 


OUoo.04 


Q SPi 

o.oo 


1 OS 
— l.Uo 


13 


1.6 10* 


33.5 


ouoo.yo 


Q fiS 
O.Do 


1 QO 
— l.OU 


13 


1.3 10* 


30.4 


PiOQA 49 
OUy4.4Z 


^ s^ 

o.oo 


114 
— 1. J.4 


J_> 


1.2 10* 


31.3 


CI no Q7 
oiuz.y ( 


l.Do 


77 
— Z. ( ( 


A 


1.3 10* 


51.3 


0110.4U 


Q SQ 

o.oo 


Qt^ 

— u.oo 


9 K 
Z.O 


7.3 10^ 


78.6 


Oo4 / .Ul 


1 fiS 
l.Do 


— 0.44 


A 


7.4 10^ 


23.3 


Dlll.Uo 


A HQ 


— u.oo 


Jj 


1.5 10* 


36.1 


CI on 1 A 

DlOU.14: 


/I 07 
4.Z / 


n Q7 
— u.y / 


R 
Jj 


2.8 10* 


22.4 


CI QQ QQ 

Dioo.yo 


/I HQ 

4.uy 


1 so 

— l.oU 


R 
13 


1.4 10* 


6.1 


CI 7p: Q7 
01 i O.O ( 


A 00 

4.uy 




— u.oo 


R 
13 


2.3 10* 


50.9 


CI 7C oo 
Dl i D.oZ 


/I OQ 


OQ 

— u.zy 


R 
15 


1.5 10* 


65.4 


fll 77 9Ci 
DIM .ZO 


1 s^ 

l.oO 


^ KA 
— 0.04 


Z.O 


4.3 10^ 


15.3 


DZU4.01 


A 00 

4.uy 


11/1 
— 1. 14 


A 

l\ 


1.8 10* 


23.2 


C07Q 9C 
DO / o.ZD 


/IIP; 

4. 10 


Sf^ 
— U.oO 


R 
13 


2.1 10* 


32.8 


CO ccn A A 
DOOU.44 


Q CO 
O.DO 


1 nn 

—i.yy 


A 

A 


9.9 10^ 


10.0 


7748.89 


3.71 


-0.39 


A 


9.5 10^ 


90.5 


7797.59 


3.90 


-0.38 


A 


1.0 10* 


80.6 


7826.77 


3.70 


-1.87 


A 


8.0 10^ 


13.1 


Y II 


log 60 = 


2.24 








5087.43 


1.08 


-0.36 


2.5 


1.310^ 


48.3 


5200.42 


0.99 


-0.65 


2.5 


1.1 10'^ 


39.9 


5402.78 


1.84 


-0.58 


2.5 


8.6 10^ 


12.8 


Zr I 


log 60 = 


2.60 








6127.48 


0.15 


-0.92 


A 


2.410" 


3.8 


6134.57 


0.00 


-1.20 


A 


2.2 10^ 


2.8 



.: NaMgAl stars 5 



5 



.5 


' 1 ' 1 


' 1 
HR448 







• 






n 
u 


• • 






-.5 


1 1 


H [ 




.5 




HR 1536 








— • **- V — • • •» • • • — 


• 




-.5 








.5 


1 1 

• 


H [ 

fi Cnc 







• 











— » 




-.5 








.5 


^ \ ^ \ 


H 1- 

20 LMi 







• • 






VJ 




• 




-.5 






— 


.5 


~: — ^ ^ ^ \ — 


H [■ 

24 LMi 




u 

-.5 


^ ^ ^ ^ 


T 

H [ 


— 


.5 





9 Com 






-.5 


» ■* - ••- «- -• 

• • 


« — 


— 




\ \ 


H ^ 




.5 




11 Aqr 


— 




• , - 







-.5 


T "••*#"• •"•B 

1 1 


_ _ _ 


— 


.5 




H [ 

HR 8472 






• 







-.5 








.5 


1 1 

• 


H [ 

51 Peg 







— ^ - - ^ ^r.#f^ T-j-n .- - \ ^ - 


-» 




-.5 


1 , 1 


1 





50 100 150 



Fig. 2. Comparison of stellar equivalent widths with those of 
EAGLNT 



in the two hottest stars and was not used for 9 Com and 
HR 8472. 

The solar equivalent widths and adopted line data are 
given in Table 2. 

3.2. Model atmosphere parameters 

Initially, the model parameters T^s, logi? and [Fe/H] were 
adopted from EAGLNT, as well as microturbulence pa- 
rameters computed from their Eq. (9). The larger set of 
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Table 2. Continued 
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observed spectral lines in this project enables the applica- 
tion of better spectroscopic statistical checks on the indi- 
vidual stellar effective temperatures, microturbulcncc pa- 
rameters and surface gravities. The initial values were then 
iteratively modified as described below. Since, for metal- 
rich stars, the temperature structures of the model atmo- 
spheres are rather sensitive to the overall metallicity, the 
model metallicities were kept consistent with the [Fe/H] 
values obtained in the analysis. 

1. Erroneous effective temperatures show up (assiiming) 
LTE as systematic dependences of derived line-by- 
line abundances as a function of line excitation en- 
ergy. To investigate such effects, the abundances of 
Fe I and Ni i lines with reduced equivalent widths less 
than log Wa (theor) / A = —5.1 (to reduce the influ- 
ence of possible uncertainties in the line broadening) 
were plotted versus the excitation energy of the lower 
level. To avoid systematic errors due to random equiv- 
alent width measurement errors, theoretical equiva- 
lent widths, VFa (theor), as suggested by Magain (1984) 
were used for this check. The effective temperature of 
each star was then modified to minimise these slopes. 

2. Similarly, inconsistent microturbulences should result 
in abundances which vary as a function of line strength 
for different lines of a species. For this test the Fei and 
Nil lines with excitation energies > 4.0 eV (to min- 
imise possible effects of errors in the excitation bal- 
ances) were plotted against log Wa (theor)/ A. The mi- 
croturbulence parameter of each star was then modi- 
fied to minimise these slopes. 

3. The surface gravities were then determined by synthe- 
sis of the spectral region near the strong Cai 6162. 17A 
line, the wings of which are sensitive to the gas pres- 
sure and therefore to surface gravity (cf. e.g. Edvards- 
son 1988). 

4. As a fourth consistency check, the ionization balances, 
which in the LTE assumption are sensitive to the ef- 
fective temperatures and surface gravities via the Saha 
equation, were studied for four elements which are rep- 
resented by both neutral and singly ionised lines in the 
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Fig. 3. Comparison of abundances derived from lines of the 
neutral and singly ionized species of the same elements, e.g., 
A[Sc/H] = [Sc/H] (lines of Sen) - [Sc/H] (lines of Sci). Stars 
classified as NaMgAl stars arc shown as filled circles, and 
non-NaMgAl stars as filled triangles. The typical effects on 
the ratios of a lowering of the effective temperature by 100 K 
are shown by the lines which end with a square symbol. Like- 
wise the effects of an increase of log g by 0.3 dex is shown by 
the crosses and of that a lowering of the model metallicity by 
0.2 dex by the stars 



investigation; Sc, Cr, Fe and Zr. While there are many 
Fei and Fell lines measured, the numbers of measured 
lines for the other species are few: only one line for 
Sci, Cri, Cm and Zrii, and at most three lines for 
Sell and Zri. Nevertheless, all the four elements show 
very similar sensitivities to the stellar parameters and 
give consistent indications of the surface gravity for 
a reasonable choice of effective temperature. This is 
demonstrated in Fig. 3. 

These four checks were then performed while the at- 
mospheric parameters were iterated until good consistency 
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was obtained for each individual star. The finally adopted 

atmospheric parameters arc given in Table 1. 

It is now of interest to investigate how these param- 
eters deviate from the values in EAGLNT, which were 
determined from Stromgrcn photometry. The new effec- 
tive temperatures are 45 K with a 1 ct dispersionof ±35 K 
hotter, which is actually somewhat less than what was 
found in an excitation-equilibrium test in Sect. 4.3.4 of 
EAGLNT. The new surface gravities are O.lTdex with a 
dispersion of ±0.13 dex higher than those of EAGLNT, 
which is within the total expected uncertainty given in 
their Sect. 4.2.1. Our [Fe/H] values differ from the photo- 
metric metallicities of EAGLNT by 0.00 ± 0.12 dex (dis- 
persion), and finally the microturbulence parameters are 
+0.05 ±0.09 km s-i higher than Eq. 9 of EAGLNT would 
give for the new values of Tes and log 5. No systematic 
dependence of these differences on the atmospheric pa- 
rameters themselves can be seen. 

3.3. Errors in the resulting abundances 

An extensive discussion of error sources (e.g., NLTE- 

cfFects and oversimplification of the model atmospheres) 
in an analysis which is very similar to the present one was 
given by EAGLNT, and will not be repeated here. We 
note, however, that we now have a much larger sample 
of spectral lines, and that the equivalent widths are more 
accurate. Furthermore, the present sample of stars spans 
a small range in metallicity and a range of only 10% in 
effective temperature, so much of the actual errors cancel 
in an analysis relative to the Sun. 

Figure 3 shows the remaining small differences be- 
tween the absolute abundances of an element as derived 
from two different states of ionization. Only for iron are 
there several lines measured from both the neutral and 
the ionized species, and the error bars show the formal er- 
rors of the mean abundances (calculated from the internal 
scatter among Fel and Fell lines, respectively) added in 
quadrature. 

Also the abundance ratios relative to Fe, [X/Fe], for 
these elements are very similar if data from neutral and 
ionized lines are used, as can be seen in Fig. 4. This con- 
firms the good quality of the equivalent width measure- 
ments. 

For C, N, O and S, the (high-excitation) lines used are 
of the neutral atom which is the dominant species due to 
their high ionization energies. When abundance ratios to 
the iron atom are plotted one finds that the ratios decrease 
systematically with increasing metallicity. If the ratios are 
instead taken with respect to the majority species Fe 11 the 
amounts of scatter around the trends diminish apprecia- 
bly. This is shown in Fig. 5. The sensitivity of the abun- 
dance ratios to effective temperature and surface gravity 
are clearly different, which, together with the smaller scat- 
ter, suggests that the ratios relative to Fell are the more 
accurate. 
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Fig. 4. Relative abundances [X/Fe] (X = Sc, Cr and Zr) 
where the open symbols represent results from lines of Xi 
and Fei and the filled symbols show ratios derived from Xii 
and Fell linos. Circles represent NaMgAl stars and triangles 
non-NaMgAl stars. The small differences between abundance 
ratios derived from neutral and ionized lines testify to the accu- 
racy of the measured equivalent widths. The effects of changes 
to Tefi and log g are indicated as in Fig. 3 



4. Results 

The resulting abundances for the nine stars are given in 
Table 3, and the line- by-line measured equivalent widths 
and abmidances are given in Table 4 (available in elec- 
tronic form). 

Most of the elements under investigation have ioniza- 
tion energies and partition functions such that the major- 
ity of the atoms are in the singly ionized state in the line- 
forming layers of these solar-type stars. When we display 
abundance ratios between different elements, we will pre- 
fer to display ratios using abundances derived from lines of 
the same ionization state in the nominator and the denom- 
inator. This procedure minimises the effects of remaining 
uncertainties in the model atmospheres, in the assump- 
tions made in the analysis and in the stellar fundamental 
parameters. The exceptions to this rule are C, N, O and 
S as discussed above. 

Fig. 6 displays our abundance results, where abun- 
dance ratios are plotted as a function of [Fe/H] derived 
from Fel lines. In the figure we have also indicated the 
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Fig. 5. Abundance ratios of C, N, O and S (derived from high excitation lines of the neutral atoms) plotted vs. iron abundance. 
The left hand panels show ratios relative to Fe i lines, and the right hand panels show ratios relative to Fe ii lines. The standard 
deviations of the abundances due to the scatter among the results for individual lines are shown by the error bars. The effects 
of lowering of the effective temperatures by 100 K are shown for one star by the open squares, and those of a doubling of the 
surface gravity by the crosses 



sensitivity of the abundance ratios to uncertainties in the 

efFective temperatures and surface gravities. The results 
for C, N, O and S were already shown in Fig. 5. 

In spite of the small number of stars the high accuracy 
in the analysis makes it worthwhile to compare the results 
with abundance trends found in other analyses. We post- 
pone, however, this discussion to Feltzing & Gustafsson 
(1997). 

Inspection of Fig. 6 shows there are no significant dif- 
ferences between the [Na/Fe], [Mg/Fe] and [Al/Fe] ra- 
tios of the NaMgAl stars and the non-NaMgAI stars. 
The six NaMgAl stars have average abundances [Na/Fe] 
= 0.06 ± 0.02 (standard deviation of the mean), [Mg/Fe] 
= -0.02 ± 0.04 and [Al/Fe] = 0.04 ± 0.01, while the 



three stars not identified as NaMgAl stars have [Na/Fe] = 
0.07±0.04, [Mg/Fe] = 0.01±0.02 and [Al/Fe] = 0.06±0.03. 
(These results use abundances, including those of Fe, from 
neutral lines.) A larger number of stars not identified 
as NaMgAl stars may be found in EAGLNT; for the 
18 non-NaMgAI stars of EAGLNT with [Fe/H] > -^0.05 
and abundances derived from EAGLNT's high resolution 
(80 000) ESO observations, we find average abundances of 
[Na/Fe] = 0.06 ± 0.02, [Mg/Fe] = 0.07 ± 0.01 and [Al/Fe] 
= 0.08 ± 0.02. Evidently, the NaMgAl stars have no sig- 
nificant enhancements of Na, Mg or Al and any possible 
enrichment of these elements is limited to ~ 0.03 dex, or 
less. 
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Fig. 6. Abundance ratios between different elements as a function of tlic iron abundance derived from Fel lines. The error bars 
show the statistical standard deviation of the mean of the abundance ratio derived from the scatter among the individual lines. 
For Sc, Cr and Zr, the plots show the mean of the abundances derived from neutral and ionized lines weighted by the available 
number of lines. For 51 Peg, the squares show the shift caused by a lowering of Tcb by 100 K, and the crosses show the effect 
of increasing the surface gravity by 0.3 dex. For Sc, Cr and Zr, these effects are shown in Fig. 4 for neutral and ionized lines 
separately. The stars are identified by symbol and [Fe/H] in the lower left corner of the figure 
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Table 3. Abundances relative to hydrogen [X/H] derived for the programme stars. A^" is the number of lines used in the analysis 
for a species. The quoted errors are only the standard deviations in the mean values due to the line-to-line scatter within the 
species with an absent error figure indicating that only one line was measured 
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In view of our failure to confirm the Na, Mg and 
Al overabundances of the NaMgAl stars claimed by 
EAGLNT, wc now briefly review the sources of difference 
between our new abundances and those of EAGLNT. Five 
primary potential sources of differences between the abun- 
dances of the two investigations may be listed: 

1. Use of additional/different lines 

2. Different equivalent widths, stellar and solar 

3. Different model parameters 

4. Different types of models 

5. Different gf values 

The two investigations use the same types of model 
atmospheres - line blanketed, flux constant, LTE mod- 
els calculated from a development of the programmes of 
Gustafsson et al. (1975) - and both use astrophysical os- 
cillator strengths derived from the lines' equivalent widths 
in the solar flux spectrum. Item four is, thus, not a fac- 
tor in the comparison, while item five is subsumed un- 
der item two. The causes of the differences between the 
two investigations are, therefore, to be found in items one, 
two and three. (A comparison of [Fe/H] abundances with 
EAGLNT shows that the present [Fe/H] values are signifl- 
cantly higher, by 0.08 dex at a mean. This is for most stars 
due to the systematically somewhat higher Tgss and log gs 
and lower microturbulence parameters of the present anal- 



ysis, and is within the range of systemic errors ascribed 
by EAGLNT to their [Fe/H] values. For 51 Peg the larger 
equivalent widths foiuid here (Fig. 2) are the main cause.) 
We now do a budget on a line-by-line basis for Na, Mg and 
Al of the sources of the abundance differences for one of 
the stars (51 Peg) which EAGLNT identified as a NaMgAl 
star. 

In Table 5 we compare our new equivalent widths and 
EAGLNT's for 51 Peg and the Sun for the Nai and All 
lines in common to our investigation and theirs; the two 
investigations do not have any Mg lines in common. A 
comparison of our new model atmosphere parameters for 
51 Peg (Table 1) and EAGLNT's reveals modest revisions, 
namely an increase of 20 K in the effective temperature, 
an increase of 0.17 in log 5, a decrease of 0.22 km s~^ in 
microturbulence and an increase of 0.11 dex in [M/H]. The 
separate effects on the Na and Al abundances of the equiv- 
alent width and log gf changes and these revisions of the 
model parameters are also shown in Table 5. For the Nal 
and All lines in common to the two analyses the new 
analysis gives [Na/H] = -1-0.25 and [Al/H] = -1-0.21, while 
for the same lines EAGLNT found [Na/H] = +0.22 and 
[Al/H] = -1-0.28. The abundance differences, for the same 
lines, between the new analysis and EAGLNT's, there- 
fore, are A[Na/H] = -^0.03 and A[A1/H] = -0.07, which 
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Line 


A [A] 


Wa [mA 


] new 


Wa [mA] EAGLNT 






A[X/H] 


(new - 


EAGLNT) 






51 Peg 


Sun 


51 Peg 


Sun 


Wa 


log g.f 




logs 




[M/H] 


Total 


Nal 


6154.23 


56.3 


39.6 


54.1 


38.2 


+0.03 


-0.02 


+0.01 


-0.01 


+0.02 


+0.01 


+0.04 


Nal 


6160.75 


74.8 


57.7 


74.9 


58.6 


0.00 


+0.01 


+0.01 


-0.01 


+0.02 


+0.01 


+0.04 


Average 
























+0.04 


All" 


8772.88 


93.8 


75.2 


109.0 


80.4 


-0.15 


+0.05 


+0.01 


-0.04 


+0.02 


+0.01 


-0.10 


All 


8773.91 


114.1 


94.4 


121.0 


98.1 


-0.06 


+0.04 


+0.01 


-0.04 


+0.02 


+0.01 


-0.02 


Average 
























-0.06 



"The new equivalent width of this line in 51 Peg is 15.2 mA smaller than EAGLNT's measurement. This poor agreement is 
caused by a nearby weaker unidentified line at 8772.54 A, which at the 60 000 resolution of the new observations is partially 
resolved from the All line and is excluded from the new equivalent width, but at the lower 30 000 resolution used by EAGLNT 
for 51 Peg was blended with the All line and was included in their equivalent width. In the Sun this nearby line is excluded 
from both the new and EAGLNT equivalent widths 



for both elements differ by only 0.01 dex from the total 
of the separate effects listed in Table 5. For the lines in 
common to the two investigations, therefore, the changes 
in the Na and Al abundances are almost entirely due to 
changes of the stellar equivalent widths, the log gf and 
the model atmosphere parameters. 

For all three Nai lines used in the new analysis the 
average [Na/H] is +0.26 ±0.01 (standard deviation of the 
mean) which is very similar to the new analysis 's value of 
+0.25 + 0.01 for the two hues in common with EAGLNT's 
analysis, while for all seven Al i lines used in the new anal- 
ysis the [Al/H] is +0.21 ± 0.01 which is the same as the 
new analysis's value of +0.21 ± 0.00 for the two lines in 
common with EAGLNT's analysis. The additional Nai 
and All lines of the new analysis thus serve to confirm 
the abundances derived from the Nai and All lines in 
common with EAGLNT's analysis. 

EAGLNT's Na, Mg and Al abundances were based on 
only two lines for each element and one of the Mg lines 
(Mgi 8712.69 A) and one of the Al lines (All 8772.87 A) 
are not of top quality so their abundances for these el- 
ements were less accurate than for elements with good 
spectral representation such as Fe or Ni. Overinterpreta- 
tion of their results for Na, Mg and Al, thus, tempted them 
to see overabundances of these three elements where none 
actually existed. The preponderance of NaMgAl stars they 
observed at McDonald Observatory compared with ESO 
is consistent with such an explanation. The 189 stars in 
EAGLNT were made up of 102 northern stars observed 
from McDonald Observatory at a resolution of 30 000, 71 

^ The 8712.69 A Mgl line is close to a stronger Pel line 
at 8713.21 A; it is not used here. The 8772.87 A All line is 
partially blended with a much weaker, but not insignificant, 
unidentified line at 8772.54 A. This All line, which at the res- 
olution (60 000) of the present investigation can be reliably 
measured in sharp-lined stars, is used here (see Table 2). 



southern stars observed from ESO at resolutions of 60 000 
and 80 000 and 16 stars observed from both observatories. 
Of the eight suggested NaMgAl stars six were observed 
from McDonald Observatory and only two were observed 
from ESO. The disproportionately large number of "de- 
tections" among their McDonald stars can be attributed 
to the larger observational scatter of the abundances, de- 
rived from the lower resolution McDonald spectra. 

Inspection of Fig. 5 and 6 shows nothing unusual about 
51 Peg compared to the other stars in these plots and so, 
although 51 Peg is metal rich, we find no sign of its com- 
position being distinctive compared to other metal rich 
stars of similar spectral type. 

Finally we briefly discuss our new determination of 
51 Peg's metallicity Our result, [Fe/H] = 0.20 ± 0.07, 
is higher than the earlier result of EAGLNT who found 
[Fe/H] = 0.06 ±0.07 and is almost identical to the results 
of Valenti (1994) and Gonzalez (1997) who found [Fe/H] 
of 0.19 ± 0.05 and 0.21 ± 0.05, respectively. Both Valenti's 
and Gonzalez's studies are also based on high-resolution 
spectroscopic observations. Most of the increase of our 
[Fe/H] compared with EAGLNT's is a consequence of our 
new equivalent widths, with smaller increases being caused 
by the higher effective temperature (AToff — +20 K) 
and lower microturbulence (A^t = —0.22) used in the 
new analysis. Valenti and Gonzalez both use very similar 
model atmosphere parameters for 51 Peg to ours; Valenti 
adopts Teff = 5724 K, log 5 = 4.30 and = 0.93 km s'^ 
and Gonzalez adopts Tcfi = 5750 K, logg = 4.40 and 
= 1.0kms-\ while we adopt (Table 1) Toff = 5775 K, 
logg = 4.35 and = 1.25kms-i. 

5. Conclusions 

Our new study does not confirm the Na, Mg and Al over- 
abundances EAGLNT claimed to find. With the aid of new 
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spectra, that for most stars have higher spectral resolution This article was processed by the author using Springer- Verlag 
and include more Na, Mg and Al lines than EAGLNT were X^T^X A&A style file L-AA version 3. 
able to use, we find that the Na, Mg and Al abundances 
of the NaMgAl stars are the same to within ~ 0.03 dex as 
those of non-NaMgAl stars of similar temperature, grav- 
ity and metallicity. The group thus appears to be illusory 
and so should be deleted from the inventory of chemically 
peculiar stars. 

The demise of the NaMgAl group of stars removes 
the basis for EAGLNT's speculation that these stars, and 
51 Peg in particular, may have stellar remnant compan- 
ions. The identification of 51 Peg's companion as a planet, 
thus, becomes more secure. We note, however, Gray's 
(1997) recent claim that the regular 4.23 d radial veloc- 
ity variation of 51 Peg is of pulsational origin, rather than 
a result of reflex orbital motion. Clearly this question must 
be settled before one can be sure that 51 Peg really has a 
planetary companion. 
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